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Abstract 

We start with a brief account of the latest analysis of the Oklo phenomenon providing the still 
most stringent constraint on time-variability of the fine-structure constant a. Comparing 
this with the recent result from the measurement of distant QSO's appears to indicate a 
non-uniform time-dependence, which we argue to be related to another recent finding of the 
accelerating universe. This view is implemented in terms of the scalar-tensor theory, applied 
specifically to the small but nonzero cosmological constant. Our detailed calculation shows 
that these two phenomena can be understood in terms of a common origin, a particular 
behavior of the scalar field, dilaton. We also sketch how this theoretical approach makes it 
appropriate to revisit non-Newtonian gravity featuring small violation of Weak Equivalence 
Principle at medium distances. 

1. Constraint from Oklo 

Around 1974 it was finally agreed that self-sustained fission reactions took place naturally in 
Oklo, Gabon, some 2 billion years ago. This is the "Oklo phenomenon." Shlyakhter |^ came to 
notice that measuring isotopic ratio of ^''^Sm left in the remnants of the "natural reactors" is 
useful to determine how much nuclear phenomena 2 billion years ago could have been different 
from what they are. He focused on the reaction n +^'*^Sm -^^^'^ Sm -I- 7, which is unique in 
that it is dominated by a resonance lying as low as 97.3 meV, nearly 7 orders of magnitude 
too small compared with MeV, a typical energy scale of nuclear physics. A very small number 
like this should be due to a nearly complete cancellation between two large effects; a repulsive 
Coulomb force proportional to a and the attractive nuclear force. 

The resonance shows up as a sharp peak in the cross section plotted against a possible 
change AE^ of the resonance energy E^. We also assume thermal equilibrium of the neutron 
flux. Even a slight change of a may result in the sizable change of E^ then of the cross section. 
Thanks to this amplification mechanism^ he derived the upper bound; |d/Q;| <,10~^^y~^, much 
better than any other results for years. 

Unfortunately, details of his derivation and the quality of the data have been left largely 
unknown. However, our recent re-analysis ||] shows that his "champion result" itself can be 
trusted reasonably well. 

To minimize the "contamination" due to the inflow from outside the reactor that occurred 
after the end of reactor activity, we relied on the latest samples collected carefully from deep 
underground. By analyzing the data on the isotopes of Sm, we computed the cross section 
(Tj^^g — (91 ± 6)kb, as bounded by the two horizontal lines in Fig. 1. Combining this with 
the shaded area corresponding to the improved estimate of the temperature (200 — 400)°C, we 
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obtained the two intersections on both sides of the peak, then the two ranges: 

{(9 ± 11) meV, Right-branch range, 
(1) 
(—97 ± 8) meV, Left-branch range. 

The right-branch range covers zero giving a "null result" as usual, while the second one implies 
that the resonance energy 2 billion years ago was smaller than today's value beyond more than 
12 a. 

We tried to see if the non-null result could be eliminated by looking at other isotopes 
^155,157qj^ii3 q^-^^ without a final conclusion yet. 

We add that the same type of analysis due to Damour and Dyson resulted in a single 
much more conservative bound as given by — 120meV < AE",. < 90meV. They could have 
separated it into the narrower two, but then with the right-branch range failing to cover the 
zero, causing a suspicion of the contaminated effect of the data. 

The bound (0) is translated into that of Aa by using AE^ — {Aa/a)Mi, with Ad^ « 
— 1.1 MeV estimated in [||, thus giving 

Aa \ -(0.8 ± 1.0) X 10"^, Right-branch range, 
a \ (0.88 ± 0.07) x 10"'^, Left-branch range. ^'^^ 

The null-result, by dividing by —2 x lO^y, gives the upper bound a/a — (0.4±0.5) x 10~^'^y~"'^, 
which "happens" to agree reasonably well with Shlyakhter's result, as was mentioned. But 
the agreement to this extent seems to be rather accidental, because, among other things, it is 
unlikely that the quality of the data available to him in the early years of the Oklo study was 
as good as ours. 

We then move on to combine our result (|l|) with those non-null result obtained recently 
from QSO, by Webb a/ Q. The error bar of the Oklo data is almost invisible on the scale 
in the upper diagram of Fig. 2a. It even appears as if a varies with time not necessarily in 
a simple manner. We argue that the time-dependence might have something to do with the 
cosmological constant problem. 

2. A small but nonzero cosmological constant and the 
scalar-tensor theory of gravity 

According to the recent detailed analyses of the Type la SN's our universe is "accelerating," 
best described by a small but nonzero positive cosmological constant, with the size expressed 
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by = A/p^j,^ ~ 0.7. The critical density p^^^ = -g^H"^ shows an over-all behavior as 
t~'^ with the coefficient which is of the order one in the Planckian unit system defined by 
c = h = Mp(= (SttG)"^/^) = 1. Note that the present age of the universe ~ 14Gy 
corresponds to ~ 10^*^ in this unit system. We then find today's value of p^^ A 10~^^°. In 

contrast, wc have a theoretical estimate A^^^ naturally ^ 1 in this unit system. The discrepancy 
of 120 orders of magnitude will be understood probably only by the "scenario of a decaying 
cosmological constant:" A is not a true constant, but decays like t^^. One of the ways to 
implement this scenario is provided by the scalar-tensor theory of gravity. 

The theory features the well-known nonminimal coupling term (l/2)-y^^(/)^i? that replaces 
the standard Einstein-Hilbert term, where 4^ = \u)\~^ shows how strongly the scalar field ^ 
couples with matter. The size of w has been best constrained by the solar-system experiments; 
tt;>3600, or ^<0.7x 10"'', though a theoretically natural estimate is of the order 1. For example, 
string theory in its simplest interpretation gives w = — 1, in contradiction with the observations. 
But the solar-system experiments may not constrain the realistic scalar-tensor theory, because 
the scalar field is likely massive, and the scalar force does not reach the distances of the solar 
system. 

Incidentally, string theory offers a nice place for the scalar field to live comfortably as a 
spinless companion of the metric tensor, and is called a "dilaton." 

A conformal transformation g^^, g^^^^ = VP{x)g^i, is crucially important in any version of 
the scalar-tensor theory, and in any theory of gravity in which a scalar field plays a role. Also 
as a terminology, under this transformation we say, we move from one conformal frame (CF) to 
another. It then follows that physics looks different from frame to frame. So the issue is how 
we can single out a correct and physical CF out of infinitely many of them. At this moment, 
we talk about only two types of them; the J(ordan) frame and the E(instein) frame. 

The J frame is the CF in which the Lagrangian takes the form with a nonminimal coupling 
term. String theory is formulated in this type of CF. For this reason, the .J frame might be called 
a string CF, or sometimes a theoretical CF. On the other hand, in the E frame the nonminimal 
coupling term has been transformed to the standard EH term. It then seems reasonable to 
start with assuming a truly constant A term included in the J frame Lagrangian, and to try to 
see what physics looks like in the E frame, expected to be a physical frame. 

Note also that in the E frame the canonical scalar field is a related to ^ by e^'^ where 
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( — 6 + Alo, and the purely constant term of A in the J frame is now converted to a term in 
the E frame; Ae'^''^'". 

We then set up the cosmological equations in the Friedmann universe with k ^ 0. We work 
out the solutions, particularly the attractor solution in the E frame. In this solution, the energy 
density p^- of the scalar field acts as A^^^, the effective cosmological constant in the E frame, and 
is shown to decay like where is the cosmic time in the E frame. This is the way we can 
implement the scenario of a decaying A. However, this smooth and monotonic behavior is not 
enough to understand the observed acceleration, which requires us to come up with something 
that behaves like a constant. To meet this requirement is not an easy task, as long as we stick 
to the idea of the decaying A. After some desperate efforts, we came to propose to introduce 
another scalar field x together with its specific interaction. Skipping the details altogether at 
this moment, we show the result 

Fig. 3a is an example of our solutions in the E frame. The horizontal axis is logt^ in the 
Planckian unit system; the present epoch is somewhere around 60. Both of the ordinary matter 
energy density (dashed curve in the lower diagram), including dark matter, and A^^ = 
(solid curve) for the energy density of the scalar fields, a and Xi fa-U off ^ t~'^ as an overall 
behavior, hence respecting the scenario of decaying A. However, we find occasional deviations 
from t^'^ , particularly the plateau behavior of p^. A nearly constant p^ mimics a cosmological 
constant, but it is remarkable that the plateau crosses p^, corresponding to fl^ — 0.5. Nearly 
in coincidence with this crossing, the scale factor (dotted curve for its logarithm in the upper 
diagram) shows acceleration, fitting the observed result. But this "mini- inflation" does not last 
forever; it dies out soon on this time scale. 

This crossing is a unique feature of the dynamics of the two interacting scalar fields, also 
indicated by the sudden changes of a and x (solid and dashed curves, respectively, in the upper 
diagram). We focus particularly on a small portion in the upper diagram, showing a sudden 
but simple change of a around the present epoch. We magnify it by the ratio 330 times in the 
vertical direction. 

The result is Fig. 3b showing a damped oscillation of a, which, like a hidden secret, is visible 
only after the huge magnification. This is a "heart" of the underlying mechanism designed for 
a small but nonzero A. The oscillation itself is so small in the amplitude that its details rarely 
affect the way of acceleration in practice. But it may show up somewhere else, particularly 
through the way a changes with time. 
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3. Varying a 



We have derived the fohowing relation between Act and Aa: 

Aa ^ a ^ , 

— =Z— CAa, 3 

where Z is a coefficient of the order 1. 

To be substituted into RHS, we take a tentatively from the above solution, and compute 
the LHS, plotted in Fig. 2b against the fractional lookback time. When compared with the 
observation, the lower diagram of Fig. 2a, the agreement appears rather poor. We have to try 
other initial values and parameters for better fits. It may take some time, though we are sure 
that we can do this without affecting much the way the universe is accelerated. In other words, 
the theory is constrained more in detail by the time- variability of a than by the acceleration of 
the universe. We still emphasize that we are coming close to understanding the two different 
kinds of phenomena in terms of a single origin, the dynamical scalar field. The most crucial 
here was (H). We outline briefiy how we reached this result j^. 

• It was necessary to modify the simplest Brans-Dicke model to be consistent with standard 
cosmology with A included, a highly non-trivial task. We also followed the scenario that 
scale invariance is broken first spontaneously so that cr is a Nambu-Goldstone boson, 
which is massless and is duly called a "dilaton." The symmetry is then broken explicitly, 
with a now as a pseudo NG boson, which is naturally massive. In this way, we have cr 
finally as non-Newtonian gravity, featuring small violation of WEP. 

• From a technical point of view we started with assuming a purely constant a in the J 
frame, then finding a cr-dependence in the E frame. The basic tool was QED. Interesting 
enough the result (^) is finite, rather than divergent, in a way of quantum anomaly, an 
intriguing concept in the relativistic quantum field theory. 

• Finally a detailed study on the choice of CF gives AG/G — — 9Z"^(Aq;/q;), and hence 
IG/G] <,10~*^^^~^^'y~^, somewhat below the current upper bounds. 
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